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ASYMPTOTIC PLANETOIDS* 



DANIEL BUCHANAN 



1. Introduction. Near the vertices of the equilateral triangles described in 
the plane of Jupiter's orbit and on the line joining the Sun and Jupiter as base 
are to be found six of the planetoids. Four of these, f Achilles (588), Patrocles 
(617), Hector (624), and Nestor (689), oscillate about one vertex in much the 
same way as the mythological luminaries of the same names circulated about the 
walls of Troy. Two planetoids are found at the other vertex, but their names 
are unknown to the author. J These six planetoids are in the vicinity of two of 
the five well known points of libration in the problem of three bodies. § The 
three other points lie on the line joining the Sun and Jupiter, one point lying be- 
tween the Sun and Jupiter, another on the side of the sun remote from Jupiter, 
and the third on the side of Jupiter remote from the Sun. One of these straight- 
line points in the case of the Sun and Earth, viz., the one on the side of the Earth 
remote from the Sun, has a physical significance in that it may account for the 
"Gegenshein."|| The equilateral triangle points of libration were considered by 
Lagrange in his celebrated prize memoir of 1772 as "pure curiosities," but re- 
cent astronomical discoveries show that these points likewise have some phys- 
ical significance attached to them. 

The object of this article is to determine orbits for these planetoids of the Sun 
and Jupiter which will approach the equilateral triangle points of libration as the 
time approaches infinity. As these orbits are asymptotic to the above men- 
tioned points of libration, we have designated the planetoids which move in such 
orbits as "asymptotic planetoids." These planetoids are considered to be of 
appreciable mass, and their perturbations upon the Sun and Jupiter are deter- 
mined. They are assumed to move in the plane of Jupiter's orbit. 



* Presented to the Society, December 31, 1919. 

t Marcolongo, II problenta dei Ire corpi, Scientia, vol. 1 (1919), No. 8. 

1 1 am indebted to Professor E. W. Brown for this information, as he communicated it 
to me when this paper was presented to the Society. 

§ Lagrange, Collected Works, vol. 6, pp. 229-324; Tisserand, Micanique CSleste, vol. 1, 
chap. 8; Moulton, Celestial Mechanics, pp. 309-318. 

|l Gylden, Sur un Cas Pariiculier du Probleme Astronomique, vol. 1; Moulton, Celestial 
Mechanics, p. 305. 
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In plain mathematical language, divorced from every astronomical applica- 
tion, the paper treats of two-dimensional asymptotic orbits near the equilateral 
triangle equilibrium points in the problem of three finite bodies. The paper con- 
cludes with numerical examples of orbits in which the ratios of the masses are 
not those of the Sun, Jupiter and the planetoids, but are chosen so that the or- 
bits near each point can be drawn to the same scale. 

Several classes of asymptotic orbits have already been obtained, but with one 
exception, number (6) below, they belong to the particular case of the problem 
of three bodies in which one body is infinitesimal and the finite bodies move in 
circles. The following is a list of the solutions which have been ob- 
tained. 

(1) Warren determined two-dimensional orbits which are asymptotic to the 
points of libration lying on the straight line joining the two finite bodies. 
(American Journal of Mathematics, vol. 38, No. 3, pp. 221-248.) 

The remaining cases were found by the author of the present paper. 

(2) Two- and three-dimensional orbits which are respectively asymptotic to 
the two- and three-dimensional periodic orbits near the straight line equilibrium 
points determined by Moulton* in the chapter on Oscillating satellites. (Amer- 
ican Journal of Mathematics, vol. 41, No. 2, pp. 79-110.) 

(3) Two-dimensional orbits which approach the equilateral triangle points of 
libration. (Transactions of the Cambrid g e Philosophical 
Society, vol. 22, No. 15, pp. 309-340.) 

(4) Three-dimensional orbits which are asymptotic to the three-dimensional 
periodic oscillations near the equilateral triangle equilibrium points which 
were determined by Buck.f This paper forms part of the paper mentioned 
in (3). 

(5) Three-dimensional orbits which are asymptotic to the isosceles triangle 
solutions^ of the problem of three bodies. (Proceedings of the Lon- 
don Mathematical Society, ser. 2, vol. 17, No. 1, pp. 54-74.) 

(6) Two-dimensional orbits which are asymptotic to the straight line equi- 
librium points when the three masses are finite. § 

2. The differential equations. Let m\, m-t, and m% represent the masses of the 
three bodies and let M denote their sum. Let a system of rectangular axes be 
chosen having the origin at the center of mass of the three bodies and the plane 
of their motion as the plane of reference. Let the axes rotate about the origin 
with the uniform angular velocity n and let the coordinates of the three bodies 
m f be denoted by (# f) y f ), i = 1, 2, 3. If the bodies are subject to their mutual 

* Periodic Orbits, chap. 5. 

t Buck, chap. 9 of Moulton's Periodic Orbits. 

t Buchanan, chap. 10 of Moulton's Periodic Orbits. 

§Rendiconti del Circolo Matematico di Palermo, vol. 45, pp. 1-25. 
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attractions according to the Newtonian law, then the differential equations 
which define their motion are 

%i — 2ny { — n 2 % { = — — -, 
Mi ox, 

y"i + 2nx' { - n 2 yi = — — (** = 1, 2, 3), 

(1) Mi oy t 

tj — U2 [ m i m 2 i w 2 w 3 m 3 m{ \ 
L »"ia r 23 r 31 J' 

r« = [fe ~ xj)* + (y { - yjY] 1/2 (i, j =1,2,3; i* j), 

where k 2 is the factor of proportionality and the accents denote derivation 
with respect to t. 

These equations admit the vis viva integral 



3 

v 



(2) £ m,<tf + y? -x?-y?) = 2U + C, 



where C is the constant of integration. 

It is shown in treatises on celestial mechanics* that the equilateral triangle 
configuration with proper initial components of velocity is a particular solution 
of the differential equations of motion (1). If the units of distance and time are 
so chosen that the mutual distances and k 2 are unity, then the angular velocity of 
rotation n must satisfy the condition 

n 2 = M. 

Since the unit of mass is so far arbitrary, it is possible, without loss of generality, 
to put M = 1. Hence n 2 is also unity. 

If the axes are initially orientated so that the #-axis is parallel to the side of 
the equilateral triangle joining m 2 and m 3 , then the values of the coordinates for 
the particular solutions are either 

Configuration I (Fig. I), 
%i = 5 (w*3 — m 2 ), y\ = - \/3 (w*2 + w 3 ), 

(3) x 2 = 5 «i + m 3 , y% = — - a/3 m h 



%3 



= — V 2 m ~^~ m )' V 3= ~ 2 ^^ mu 



* See Moulton's Introduction to Celestial Mechanics, pp. 309-11. 
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(4) 



*i = 2 ( ms ~ m ^< 



%=-„»i + m z , 
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Configuration II (Fig. 2), 
1 



y\ = — s V3 (w 2 + m 3 ), 
J2 = g V3 mi, 



*3 = - f g w«i + ^2 J, ^a = 2 V3 mi. 



[June 




Fig. 1. — Configuration I. 



Fks. 2. — Configuration II. 



Since the origin of coordinates is at the center of mass, it follows that 

m\%\ + vn^H + m&,% = 0, m^y\ + m^yi + m»ys = 0, 

and therefore the coordinates of one of the bodies may be eliminated from the 
differential equations. Let us suppose x 3 and y 3 are thus eliminated by 



(5) 



= {m\X\ + nhx^), 



yz- (w^i + tmy*). 

m 3 



When these substitutions are made in (1), and k 2 and n are both put equal to 
unity, the differential equations of motion become 
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x 1 — 2y x — Xi = 



m 2 (xi — Xi) (wi + nt3)Xi + m 2 x 2 

~3 



Pl2 Pl3 



7l 1" &X X Ji — — — , 



x 2 — 2y 2 — x 2 = 



„ 3 .3 

Pl2 Pl3 

_ _ m\(xj — aci) _ m\X\ + (m 2 + m 3 ):*; 2 



Pl2 P23 



(6) 



y 2 -t- z# 2 — J2 — — 5 — 5 , 



P12 P23 



P12 = [(*i - * 2 ) 2 + (yi - y2) 2 ] 1/2 , 

= — [{(1 — «* 2 )*i + ^2^2 } 2 + {(1 — m 2 )yi + m 2 y 2 } 2 } 1/2 , 



P13 

w 3 



— .[{niiXi + (1 — mi)x 2 } 2 + [ntiyi + (1 — Wi)y 2 } 2 ] 1/2 



wj 3 



3. The equations of variations. Let 

(7) 

^2 = „ mi + w 3 + eu 2 , y 2 = — - s/Z mi + ev 2 , 



where U\,. . ., v 2 are new dependent variables, and e is an arbitrary parameter. 

The additive quantities eui eu 2 denote the components of displacement of m\ 

and m 2 from the vertices of the equilateral triangle in configuration I (Fig. 1). 
We shall not consider configuration II (Fig. 2) in detail as the coordinates in 
Fig. 2 differ from the corresponding ones in Fig. 1 only in the sign of \/3. Since 
\/3 does not occur explicitly in the differential equations, we may change the sign 
of \/3 in any solutions which we obtain for configuration I and thereby obtain 
the corresponding solutions for configuration II. 

When equations (7) are substituted in (6) the right members of the differen- 
tial equations can be expanded as power series in e which will converge for e 
sufficiently small. Since the expressions in (7) that are independent of « are 
particular solutions of equations (6), there will be no terms independent of e after 
the substitutions (7) have been made. Hence the factor e can be divided out of 
the equations, and when m$ is replaced by 1 — mi — m 2 the differential equa- 
tions (6) become 
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(d, - fjm - ( 2D + ^ (1 - 2mA Vl + 0u 2 - ^ m 2 v 2 

= eUf + e 2 U[ 3) + • • • + e"- 1 ^ + . • ., 

i 2D - ^p (1 - 2m 2 )\u 1 + (d 2 - |) Vl - ^0 m 2 u 2 + 0v 2 
(g) = V3 t€V?> + eW[ 3) +..- + e n - l V™ +■■■], 

9 3\/3 , / n . „,9 \ / on 3a/3 \ 

— - WiMi — —V- «tii»i + I D 2 — 3 + - mi ) m 2 — I 2D — — — wi I v 2 
44 \ 4 / \ 4/ 

= «C7f + « 2 C7f + • . - + e n - 1 Ui n) + ■■■, 

3\/3 .9 ,L. 3\/3 \ , / n2 9 \ 

— — t— w^Mi + - ntiVi + [2D H — mi] u 2 + [D 2 — - mi] v 2 

4 4 \ 4/ \ 4/ 

= V3 [eVf + £ »V<» + • • • + e"" 1 ^ +•••], 



where D denotes the operator d/dt, and t/f',. . ..V^ (k = 2 «,...) are 

homogeneous polynomials in Mi, %, »i, % of degree k. In all these poly- 
nomials, V^3 occurs only as a factor of the odd powers of vi and v 2 , considered 
together. 

If we consider only the linear terms in (8) we obtain the equations of variation. 
They are 

(d 2 - fj «x - 1 2D + 3 -^5 (1 - 2m 2 )\v i + 0% - 3 -^2 = 0, 

| 2 D - ^p (1 - 2m 2 )l % + (d 2 — |) ^i — ^ »»*«» + 0% = 0, 

i ^— wii^i + f D 2 — 3 + 2 w i ) M 2 — ( 2D — - • mijv 2 = 0, 

3\/3 ,9 L L, 3\/3 \ , / n2 9 \ 
j- wiMi + - mi7)i + I 2D H — mi)u 2 +iD 2 —-mi jv 2 = 0. 

4. The solutions of the equations of variations. The solutions of the equa- 
tions of variation can usually be obtained by differentiating the generating solu- 
tions, equations (3), with respect to the arbitrary constants which denote the 
initial time and the scale factor.* This is not possible in the case under con- 
sideration since the generating solutions are constants. 

Equations (9) are linear differential equations with constant coefficients, 
and in order to obtain their solutions we consider the operator D as an alge- 



W 9 



* Poincarl, Les Mithodes Nouvelles de la MScanique CSleste, vol. 1, chap. 4. 
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braic quantity and equate to zero the determinant formed from the coefficients 
of «i, u%, vi, i> 2 in (9). This gives 



(10) 
A = 



D*-\, - 2D - ^ (1 - 2m,), 0, - ^ m, 
2D _ 3V3 (1 _ 2w0> D2 _9 _3^/3 W2; Q 

9 3\/3 r>2 q . 9 on . 3 V3 

- 7 «i, — — — Wi, A? 2 — 3 + r Wi, — 2L> H — «i 

4 4 4 4 

3V3 9 on , 3\/3 n2 9 . 

— — - — mi, - mi, 2L> + — — - mi, D l — - m\ 

4 4 4 4 

This determinant reduces to 



= 0. 



A = D- 



(11) 



D 6 + 2£> 4 + < 1 + — (W1W2 + w 2 w 3 4- W3W1) f -D 2 

27 1 

+ t (W1W2 + w 2 w 3 + W3W1) = 0. 



If we neglect the factor D 2 , the preceding equation is a cubic in D 2 and its dis- 
criminant is 

27 

— (wiw 2 + mim 3 + W3W1) 2 [27(wiw 2 4- wimi + m 3 mi) — 1]. 

Since Mi + fw 2 + W3 = 1 and since mi, w^, and w 3 are all positive, this dis- 
criminant is positive, and therefore the cubic has one negative root and two con- 
jugate complex roots. Let these roots be X 1( X 2 , and X 3 where 

(12) Xi = V 17 ! *> X 2 = m + V^l v, X 3 = ». - V^l v. 
Then the solutions of (11) are 

(13) D = 0, 0, ± Xi, ± X 2 , ± X 3 , 
and the solutions of the differential equations (9) are therefore 

mi = on + a 2 t + a 3 e M + a ie ~ M + a^" + me'™ + we™ + a 8 e" X3 ', 
vi = 0<% + 0<°W + /3< V - 0Pa#-™ + 0<V ( - ^^e^' 

(14) % = Tl 0) «i + 7 2 °W + TiV + yi l) a ie - M + 7 i V + y^ase-™ 

+ 7 ?W s( + yfase-™, 

v 2 = «i 0) oi + 5 2 °W + S^/' - B?a#- M + 5< V - 5£ 2 W -X! ' 
+ Specie™ - 8?W -X *. 
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where en, . . . , as are the constants of integration. The remaining symbols /3j , 
. . . , sf denote constants which are defined by the following equations and prop- 
erties : 

pf > : — f- - ^ (1 - »m) - - (m. - v 7 3 m 2 ) 

1 — mi [9 3 3 

+ \wi\{mi +W3) + 4whmz\ I, 

9(wi - 1) J 



(0) _ ( w i — 1) ("\/3 wi + 9miw 2 ) + \/3 \m\{wii + nt 3 ) + 4w 2 w 3 } 

9m 2 (wi — l) 2 



«<<» = - — r 

2m 2 L 



Q (a) 4. g^ 3 - W2) \/3" 

9(W! - 1) 3 



(1 - 2m 2 ) M 0) + ^ ^ + ^ , 



,(0) _ m 3 — m 2 ^) _ wi *(o) _ mi + 2w 3 

— ,_ j / 2 — J 2 / ' 

V3(>Mi — 1) mi — 1 \/3(wi — 1) 



J 2 



< = 3\/3 [(m 3 - m 2 )Xf + 2\/3(l - wtQXj + 3(m 2 - m 3 ) ] 
4Xj - Xf (5 - 9mi) + Qy/3(m - m 2 )X,- + 27(1 - mtf 



T « = 2V3 A f _ 9\ ^ + 4V3 ^ _ l_-a* 
9m 2 \ 4/ 9m 2 2w 2 

5 <0 = 2 V| (xf - 3 -) - { 4 -^? X,. + L^-^?} fi> (i = 1, 2, 3). 

9m 2 \ 4/ v 9m 2 2m 2 ; 

The following pairs of constants differ only in the sign of \/3 : 

0i°, z^ ; 7i°,7i°; Si 4) .#; (* = 1.2,3). 

The following pairs are conjugate complex : 

Pi » ~ P2 > 7l . 72 . "l > — °2 . Pk > Pk , Ik » 7* , °k > °k > \K — i. ^; 

The solutions for D in equations (13) are called the characteristic exponents. 
Since the original differential equations (1) admit the integral (2), then two of 
the characteristic exponents will be zero.* It was to be expected, therefore, 
that two roots of (11) should be zero. Poincaref has also shown that if equa- 

* Poincarfe, loc. cit., p. 188. 
f Loc. cit., p. 69. 
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tions (1) do not contain t explicitly, which is always the case in all problems of 
mechanics in which there is a conservative system, then the characteristic expo- 
nents are always equal in pairs but opposite in sign. It is this property of the 
solutions that makes the determinant (11) even in D. 

5. Construction of Asymptotic Solutions in e~ X2 ' and e _Xs '. We shall now 
construct solutions of the differential equations (6) which are asymptotic in the 
sense of Poineare, * that is, each term of the solution must have the form 

e" P{i), 

where c is a constant having its real part different from zero, and P(t) is a periodic 
function of t or, in particular, a constant. Such solutions will therefore approach 
zero as t approaches + oo or — oo according as the real part of c is negative or 
positive, respectively. 

The only terms of the solutions of the equations of variation which are 
asymptotic in the sense just defined are those in e* ! ' and e* a . These terms 
approach zero as t approaches — oc or + oo according as the + or — signs are 
taken with X2 and X3. 

In this section we shall construct the solutions which approach zero as t ap- 
proaches + 00. These solutions will obviously involve powers of e~ !/ and 
e~ . In the next section we shall show how the solutions in e '' and e *' 
can be obtained directly from the solutions in e~ 2( and e~ a . 

Only the formal construction of the solutions is considered in §§ 5 and 6, but 
the convergence of the solutions obtained is established in § 7. 

We propose to integrate equations (8) as power series in « and this is the reason 
why e was introduced in equations (7) . Accordingly we substitute 

(15) u { = ^ Uij e j , »,■ = ]£ v {j e j (i = 1, 2), 

} = j = 

in equations (8), and as we shall have frequent occasion to refer to the resulting 
equations we shall cite them as (8')- These equations (8') are to be satisfied 
identically in e and we may therefore equate the coefficients of the same powers 
of e. In this way we obtain sequences of differential equations in u ti and v {j 
which can be integrated step by step, as we shall show, subject to suitable initial 
conditions and restrictions. 

The first two steps of the integration will be considered in detail, and then an 
induction to the general term will be made to show that the process of integra- 
tion can be carried on indefinitely. 



* t ™ cit, p. 340. 
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Step 1. Terms in (8') independent of e. The terms in (8') which are indepen- 
dent of e are obviously the same as the equations of variation (9) if the subscripts 
on the dependent variables in (9) are altered so as to read «io, vw, M20 and % 
instead of u u u- 2 , vi, v%. Hence the solutions of these equations are the same as 
(14) with the corresponding changes of subscripts. Now only two of the ex- 
ponentials, viz., e~* 2 ' and e~^, approach zero as t approaches + 00 and we 
therefore put equal to zero the arbitrary constants associated with the other 
exponentials. The solutions at this step are, then, 



_ .„«>) „-™ J. „(0) „-M „.„ = _ „(2) „(0) a -\*t _ ft (3) „(0) ~Ml 

(16) 



w 10 = a?' e "" + <' e "", D10 = - PI" <' e ™ - ft" «« ' e 



«*» = 7?> «f e"" + 7^ af <T X ", v w = - 5< 2 > a<°> «"" - 5< 3 > a™ e~™ 

where «g 0) and a^ are arbitrary, and p^ 21 , . . . , df^ are defined as in equations (14). 

It is evident that the complementary functions at all the succeeding steps of 
the integrations will be the same as (14) and that after the exponentials are re- 
jected which do not have the proper form, two constants of integration will re- 
main undetermined. Hence, besides rejecting all the exponentials except e~ 
and e~ ut , it is necessary to impose two initial conditions upon these solutions 
so as to determine the arbitrary constants arising at each step of the integration. 

Let us suppose that 

(17) «i (0) = e lt uz (0) = 7, 

where «i and 7 are arbitrary. Since ui and m 2 are multiplied by the arbitrary 
parameter « in (7) we may put t x or 7 equal to unity. Let us suppose «i = 1. 
When the conditions (17) are imposed on (15) it follows that 

(18) «io (0) = 1, My (0) = 0, 

«20 (0) = 7, % (0) = (j = 1 00). 

Now applying these conditions to the solutions (16) we obtain 

j« _ v™ - y „(o) = 1 - ^ 2) 

6 y?-y?' ' 7^-7f 

Since y™ and y® are conjugate imaginaries, it follows that «g 0) and a® are like- 
wise conjugates. 

In order to unify the notation, we put 

(0) _ „(10) _ a(t) .,(0) _ o(10) „(2) (0) _ (10) _ .5(2) (0) _ o(10l 

flOA 6 _ 10 ' P2 a 6 — PlO > 72 a 6 a 20 . °2 a S P20 1 

^ „(0) _ „(01) _ 0(8) (0) _ o(01) (3) (0) _ (01) _ g(3) (0) _ fl (01) 

«8 a 10 > Pi a i ~ - PlO 1 72 a 8 a 20 > °2 a S P20 • 
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Since a ( «\ a^; $\ /3 2 3) ; y?\ yf; 5 2 2) > S 2 3) are conjugate pairs, it follows 
that the symbols in (19) which differ only in a permutation of superscripts are 
likewise conjugate pairs. This notation will be adopted in the sequel to 
denote conjugate complexes. Thus a tl and a,[ (;' 9^ k), are conjugate 
imaginaries. 

On employing the above notation, we obtain for the desired solutions at this 
step 

(20) u i0 = «}> 0) e~™ + aft" e- x ", v i0 = /#•> e~ M + &» e~™ (i = 1, 2), 

where a;J 0) , . . . , /3JJ 1 ' are linear in 7. 
If we put 

x 2 = ii + V— 1 ". X3 = m — V— 1 "» 

and suppose that 

4 10) = i la$ + V- 1 &SM, 0.r = I [4° + V~l $>]. 

then the solutions (20) become 

u i0 = e~"' [all ) cos vt + fc,^ sin vt], 
v i0 = <r"' [c$ cos v* + dj? sin vt]. 

This second form of the solutions is more convenient than (20) for numerical 
computation but is more cumbersome in obtaining the solutions at the succeed- 
ing steps. 

Step 2. Coefficients of e to the first degree. When the solutions (20) have been 
substituted in (8')» the differential equations obtained by equating the coeffi- 
cients of e to the first degree are 

(d* - I) an - Ud + 3 -^? (1 - 2m0 } v u + u 2l - ?^§ miVil = U 1U 

Ud - ^ (1 - 2wt 2 ) I «„ + (d 2 - I) i* - ^J ***, + vn = V11, 

( 22 ) - I WiMi, - ^ «n, + ( I? 2 - 3 + I mi) «* 

- ^2D - ^ Wl ) *, = J7 n , 

~- WiMn + - wi^ii + (2D + -^— wi J W21 + f D 2 — j m A "21 = ^»: 
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I/u = A[\ 0) e~ 2 ™ + A\\ l) e~^ + Xs)( + 4<f <T 2Xi ', 
V n = B™ e~ 2M + B™ e~ (M + Xs)( + £<?> <r 2X3/ , 
I/„ = A ( 2 \ 0) e- 2X « + A%> e~^ + x ' )( + A 2 f e~ 2U , 
Vn = B' 2 \ 0) e- 2 * 1 + B™ e~^ + X3) < + B<»> e~ 2 ^, 



where A™, •••> B^i are quadratic expressions in 7 having constant coefficients 
The constants like A™ and A™ which have the same subscripts but which have 
their superscripts reversed are conjugate imaginaries. The other constants 
A[\ v ', ...,B<i° are real. 

The complementary functions of equations (22) are the same as (14), and the 
particular integrals can be found by the method of the variation of parameters 
as at the previous step. These particular integrals are 



(24) 



where 



(25) 



«*i = £ <#* e ~^ + * X3) ' (J + k = 2, i = 1, 2), 

i, k = 



3, k = 



(jk) e -<.fa + k\ 3 )t 



A 0'*) A 0'« 

Jjk) _ _^1 o(jk) _ &2 



11 "" A[- (j\ 2 + k\ 3 )]' ' n A[- (/A* + feX,)]' 

x (Jk) _ ^3 aO - *) = ^4 

21 A[- (/X 2 + feX,)]' 21 A[- (;X 2 + feX,)]' 



The preceding symbols are defined as follows : 

A[ — (/X 2 + k\s)] denotes the determinant A in (10) when D has been re 
placed by — (/X 2 + &X 3 ) ; 

^i k) , i — 1» 2, 3, 4, denotes the preceding determinant A[ — (;'X 2 + k\ 3 )] 
when the elements of the ith row, reading from top to bottom, have been replaced 
by A™, B™, A 2 f , B%\ respectively. 

The constants in (25), viz., a u , . . . , |3 2 i are found to have the same form 
as the constants in (23). Thus a u and a^ are conjugate complex if / ^ k, 
but real if / = k. 

If the right members in (22) contained any constant terms, or terms in t or in 
e ±x *', i = 1, 2, 3, that is, terms which are the same functions of t as the terms in 
the complementary functions, then the denominators in (25) would vanish and 
the particular integrals would not take the form (24), but would contain terms 
in t, t 2 or t £ ±x »'. Hence, such terms would appear in the particular integrals 
if and only if the right members contained constants, terms in t or e ±x ' ( , re- 
spectively. In order, then, to show that there are no terms in the particular 
integrals which do not satisfy Poincard's definition of asymptotic solutions- it 
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is sufficient to show that the right members do not contain any terms which 
are the same functions of t as any part of the complementary functions. 

To obtain the complete solutions of (22), we combine the complementary 
functions, equations similar to (14), and the particular integrals (24). The par- 
ticular integrals have the desired form for asymptotic solutions, but we must re- 
ject all the terms of the complementary functions except those in e~ M and 
e -x»( ^j s j g p 0Ss jtj| e> f course, by equating to zero all the arbitrary constants 
of integration except those associated with the exponentials e~ M and e~ Xs/ . The 
desired solutions are therefore 

Uu = a W e ~M + a (» ,-W + a J10) e ~2M + a (ll) ,-©. + X.X + a (02) g -2X 3l> 

va = -pf ' «<" e~™ - pf > «<» e'™ + fl'f e~ 2M + 0™ e~' M + x ° )( 

+ /3'f <T 2Xj( , 

u*. = y? «?' e~™ + 7?' «s U e~ Mt + ag 0) e~ 2X *' + ajj 1 ' <T (X * + Xs j ( a<f ^ 2Xsi , 

where a s and a^ 1 are the constants of integration. By virtue of the initial 
conditions (18), these constants must satisfy the equations 

«?> + <4° = " kf + ail "+ "if] - ^\ 

or 

(1) _ *<*> 7 <«> - g^ (1) = *<» - g^ T<*> 

7 W- 7 » ' *° 7 < 3 >- 7 < 2 > * 

Since a^, a[\ 2) and af®, a[\ 2) are conjugate pairs and a[\ l) and a^ 1 ' are real, 
R^ and i?^ are likewise real. But as 7 £ s) and y^ are conjugates, it follows 
that aj and «g are also conjugates. 

In order to unify the notation, as at step 1, we put 

„(1) _ „(10) _ o(2) (1) _ fl(lO) (2) (1) _ (10) _ g(2) (1) = fl (10) 

a 6 — a ll ) P« a « ~~ Pll > 12 "6 "21 > °2 "6 M21 > 

„(D _ ,,(01) _ fl(3) (1) = fl (0l) (3) (1) _ Ull) _ g(2) (1) _ fl (0l) 

#8 — a ll Pi a 8 Pll > T2 "8 "21 > °2 "8 P21 > 



and the desired solutions at this step take the form 

(26) lk = 

Vn = J2 P™ e ~° M + * X,) ' (j +k = 1 or 2, * = 1, 2). 

j, k = o 
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If we put 

^2 = M + "\/ — 1 v, X 3 = n — v"— 1 y, 

and if we suppose that 

(27) J ^ 

^ = \ [4 k) + v=i d^i $» = \ [c^ - v~i <&*>] a = i, 2), 

where b^ ) and d fl are different from zero when / =^ k, but equal to zero when j 
= k, then the solutions (26) become 

u h = e-"' [a£ 0) cosvt + bjl 0) sin vt] + e~ 2lA Q og^ + <$ 0) cos 2v * 

+ b% 0) sin 2» t\ (« = 1, 2), 

"a = e~ M< t4; o) cos » * + 4 10) sin r t] + e~ 2 >" Q 4 n) + c^ cos 2v * 

+ 4 20) sin 2v t\. 

The succeeding steps of the integration are similar to the preceding step, and 
an induction to the general term will now be made to show that the process of 
integration may be carried on for any desired number of steps. 

Let us suppose that u u , v it , i = 1, 2, have been computed for / = 0, . . ., 
n — 1, and that 



i + i 



«« = J2 a $ 



(»*) --&>« + **»)« 



e 



j,k = 



(28) 

"« = 'l>f } ^ 0> + * X3)< 0' + * = 1, 2 Z + 1). 

j, k = o 

where a^ and a,- ; ; /?# and $/ are conjugate complex if j^kbut real if ;' = k. 
In order to make the induction it will be necessary to show from the differential 
equations in u in and v in obtained from (8') that the solutions for these variable 
are the same as (28) if I = n. 

Step n + 1 : Coefficients of e" in (8') • Let us consider the differential equations 
obtained by equating the coefficients of e* in (8') after the various solutions in 
(28) have been substituted. They are found to have the same form as (22) if the 
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second subscript on the variables and the right members is changed to n. The 
right members, however, have the form 



(29) 



n + 1 

TT - V AW 0-^* + * Xs)i 
*« / j *" ' 

j. k = 

V fa = J2Bg k) e-<* + m)t (j + k = 2, 3 n + 1; * = 1, 2), 

j, k = o 



where .4^ and 5^ are complex quantities similar to those in (28). The com- 
plementary functions of these differential equations which have the desired ex- 
ponentials are the same as (16) if the second subscript on the variables and the 
superscript on the constants of integration are replaced by n. The particular 
integrals of these differential equations can be found as at the previous step by 
the method of the variation of parameters. They are 



= y eta ,-<* + »*< s y a w» e - 



_ y-* A 2« *>~°' X2 + kX3)t = V* flW fe) e _0 ' X! + ftXs >' 



(30) n + l A 0' + t) « + l 

_ y^ A 3« e -oxi + *x.)« = y^ 



w _ y» "3« -yx 2 + *x,x ^ y-» uk) -0x2 + *X!)< 

2W . A^. A r _ cv\. a. fc\.M ,4-J* n > 



« + 1 A y + *) « + 1 

„ _ 'y 1 A 4n t e -yx ; + ku)t _ y-» a<Jk) e -U\* + k\ 3 )i 

2 " ,- £ A[ - 0'X 2 + feX,) ] ,• j^ 

0' + k = 2, 3, ...,n+ 1). 



The various A's in the preceding equations are defined as follows : 

A[ — (;X 2 + £X 3 )] denotes the determinant A in (10) if D is re- 
placed by — (/X 2 + &X 3 ). Since / 4- k is not less than 2, these de- 
nominators do not become A( — X 2 ) or A( — X 3 ) and are therefore different 
from zero. 

Ai„ \ i = 1, 2, 3, 4, denote the preceding determinant if the elements of 
the ith row, reading from top to bottom, are replaced by A^„ , B^ n , A^ n and 
B 2n , respectively. 

On determining the constants of integration .a" and ag by the initial condi- 
tions (18), and then unifying the notation as at step 2, the complete solutions 
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which have the desired form are found to be the same as (28) if I = n. This 
completes the induction. 

The solutions, therefore, at the general step n are 

i, k - 

(31) -V#'<.-^ + JXl " 

"in ~ / j Pin " 

(i = 1, 2; ; + ft = 1, 2, . . ., n + 1; n = 0, 1, 2, . . .). 

When equations (31) are substituted in (15), the asymptotic solutions of the 
differential equations (8) are then found to be 

(32) «-ii.*-o 

00 « 

€ "« = £ £ ^ e ~ ux ' + m)t e " (* = L 2 ; / + * = 1. 2, . . .,«), 

« = 1 j, ft = o 

where the second subscript on the constants af® and /3,^ has been made to 
conform with the powers of e. 

If these solutions are expressed in trigonometric form by the substitution 
used in (27) we obtain 

« n + 1 

««*=/< (+ vs. - o - £ 2 E' £* [{«••«' ~ h2k ~ *' c ° s (2 * - d ** 

/oo\ n = 1 J = 1 fc = or 1 

+ 6g* - *■ 2 * - 1J sin(2£ -l)rt) e~ (2j ~ *>" + [a$- 2k) cos 2k v t 

+ 6,^ 2 * ) sin2fe^} e- 2j >"]e n , 

««j = ft (+ V3, - (*'» = 1 ' 2 )' 

where 2' denotes that the highest value of 2Z — 1 or 2Z is n, and 2* denotes that 
the lowest values of 2k — 1 and 2k are 1 and 0, respectively. The function 
gi ( + a/3, — t) is the same as /,• ( + a/3, — if the constants a and 6 with their 
subscripts and superscripts are replaced by c and d, respectively, with the same 
subscripts and superscripts. When these equations are substituted in (7) we 
obtain the parametric equations of the orbits which approach the vertices of the 
equilateral triangles in configuration I, viz., 



1922] ASYMPTOTIC PLANETOIDS 425 

xi = - (ra 3 - w 2 ) + A ( + a/3, - t), 

•\/3 /- 

^i = — (nt 2 + m s ) + gi ( + V3, - t), 



(34) 



ft = 5«i + »(j+/!( + V3, — t), 

\/3 A" 

y 2 = — -^- wi + g 2 ( + V3. - 0. 

x 3 = (wi *i + w 2 x 2 ), 

ys = (mi y x 4- w 2 ;y 2 ). 

m 3 

The orbits which approach the configuration II are obtained by changing the 
sign of a/3 in (34). They have the equations 

xi = - (w s - w 2 ) +/i (- V3, — 0. 

a/3 
Ji = - ^r- (wg + wa) + gi (- V3, - «), 



(35) 



#2 = g m i +" w s + h (~ V3, — 0. 

y 2 = ^- wi 4- g 2 (- v'3, - t), 

%3 = — _ (tw, ^ -f. W2 ^ 2 ) 
m 3 

j3 = (wi yi + m 2 j 2 ). 

*M 3 

Equations (34) and (35) therefore represent the orbits which approach the 
vertices of the equilateral triangles in configurations I and II, respectively, as 
the time approaches plus infinity. These orbits deal with the future of the system 
and, at the risk of being censored, we shall cite them as future I anidfuture II, 
respectively. 

6. The asymptotic solutions in e M and e Mt . This section deals with the 
asymptotic orbits which approach the vertices of the equilateral triangles as 
the time approaches — oo . Such orbits deal with the past while the previous 
orbits forecast the. future. So then the mathematical astronomer can say with 
more than poetic license 

"Backward, turn backward, O time in your flight." 
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It is obvious that these orbits can be constructed by making the same use of 
the exponentials e '' and e ** as was made of e~ * and e~ a in constructing 
the previous orbits. This construction will not be considered in detail as we 
shall show that the past orbits can be obtained from the future orbits in a very- 
simple way. This method, it will be observed, is the converse of the adage of 
"history repeating itself." 

Let us consider the differential equations (8) and let their solutions (33) be 
denoted by /<( + y/z, — t), g t ( + y/z, — t),i = 1, 2. Next consider the effect 
of changing the signs of t, y/% v t in (8). The left members of the first and third 
equations of (8) remain unchanged while those of the second and fourth change 
signs. Since u[ and v[ are even in \^3, »i, and %, considered together, the above 
changes of signs will leave these expressions unaltered. But as the right members 
of the second and fourth equations in (8) contain the factor v3, these changes of 
signs will produce a change of sign not only in the left members but also in the 
right members of these two equations and the minus sign can be cancelled off 
in both equations. Hence the differential equations (8) are unchanged if the 
signs of /, y/% Vi and % are changed. These changes of signs have no effect upon 
the initial conditions (18) and consequently if we make the same changes of signs 
in the solutions /,( + \^3, — t), g,( + y/Z, — t),it will still leave them solutions 
of (8). Thus/<( - s/Z, + t) and - g,( - \/3, + t) are solutions of (8) and 
the corresponding solutions of (1) are 

a/3 /— 

y\= ~ (nh + m 3 ) — gi (- V3, + t), 

od = 5 mi + m 3 + / 2 (— V3, + t), 

m V3 

y* = — ^ mi - gi (- s/3, + t), 

x 3 = (mi Xi + m 2 X2), 

m 3 

yz- (mi yi + m 2 y 2 ). 

These are the orbits which approach the vertices of the equilateral triangle in 
configuration I as the time approaches — 00 , and will be cited as past I. The 
corresponding orbits for configuration II, past II, are obtained from (36) by chang- 
ing the sign of y/z. They are 



1922] ASYMPTOTIC PLANETOIDS 427 

%i = o ( w s> - m *) + /i (+ V3. + Oi 



?i 



= - ^ (mt + m 3 ) - g, (+ V3. + 0. 



(37) 



1 

*2 = 5 wi + m 3 + / 2 (+ V3, + 0i 

x 3 = (mi Xi + tfh X2), 

Wit 

y 3 = — — (wi yi + m 2 yi). 
m 3 

From the form of the equations of these four orbits' it is evident that the past 
II orbits are obtained from the future I orbits by changing the sign of t in the 
latter and reflecting in the x-axis. The same relation exists between past I and 
future II orbits. 

7. The convergence of solutions. Only the formal construction of the solu- 
tions (8) has been made, and we shall now consider their convergence. 

The convergence of the solutions (32) of the differential equations (8) de- 
pends upon the form of the characteristic exponents in terms of which the so- 
lutions are expanded. Now it has been shown by Poincare* that such solutions 
as (32), (34) and (35) will converge as t approaches + oo provided that the real 
parts of the characteristic exponents in terms of which the solutions have been 
expanded are different from zero and positive, and likewise solutions such as 
(36) and (37) will converge as t approaches — <x> provided that the corresponding 
exponents are different from zero and negative. The characteristic exponents 
of this problem are 0, 0, =±= Xi, =±= X2, =*= X 3 , equations (13), but only X 2 and X 3 
have their real parts different from zero, and it is in terms of only these exponents 
that the solutions have been expanded. Hence the solutions (34) and (35) 
converge as t approaches + 00 , and the solutions (36) and (37) converge as / 
approaches — 00 . 

8. Illustrative examples. We shall conclude this paper with illustrative 
numerical examples. The values of the masses chosen are not the ratios of 
the masses of the Sun, Jupiter and the planetoids mentioned in § 1, but more sim- 
ple values to illustrate the nature of the orbits. 

Let mi = 0.2, w 2 = 0.3 and m% = 0.5. Then M, the sum of the masses, is 
unity, and n, the mean angular motion, is + 1 or— 1. We have chosen the 
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counter-clockwise direction of rotation and put n = + 1. The vertices of 
the equilateral triangle for configuration I are (0.1, 0.4\/3), (0.6, — 0.1\/3) 
and ( - 0.4, - 0.1\/3), and for configuration II (0.1, - 0.4\/3), (0.6, 0.1\/3) 
and ( — 0.4, 0.1\/3) for mi, m* and mz, respectively. 

The computation has been carried out for only the linear terms in t for the 
orbits (32). The values of the various constants together with the equations 
which define them are listed in the following table. The conjugates of the var- 
ious terms are omitted. 

Tabls I 



Constant 


Equation 


Value 


A 


(11) 


252(£>e+2Z?<+3.0925Z? 2 +2.0925) 


Xi 


(12) 


V-. t 


Xj = /» + V— i " 


(12) 


0.688+0.987 V-l 


4 2) 


(14) et seq. 


0.152+0.948 V-l 


y? 


(14) et seq. 


-0.246+0.547 V^l 


«, w 


(14) et seq. 


-0.630-0.315 V-i 


4 0) = -8? 

_ g(2) (0) = g(l0) 
D 2 "6 °10 


(16), (19) 
(16), (19) 


0.500-(0. 225+0. 914-y) V^l 
-(0. 289+0. 866y)-(0. 440-0. 139v) V^l 


v («) a <0) = a d0) 
T 2 6 20 
_j(2) (0) = fl (l0) 
°2 a 6 "20 


(16), (19) 
(16), (19) 


0.500v+(0. 328+0. 225r) V-l 

(0. 386+0. 288-j-) + (0. 015-0. 578t) V~-l 



The solutions for t u u t vi, t ut, and e ih are 

««i = «-0.688/ [ cos 0.987* - (0.450 + 1.828-y) sin 0.987*] <H , 

61 = e-0.688t [(0.578 + 1.732-y) cos 0.987* + (0.880 - 0.278 7 ) sin 0.987/] « + ••-, 
««2 = e-0.688< [ y cos 0.987* + (0.656 + 0.450-y) sin 0.987*] « + ••-, 

62 = e -0.688( [(0.772 + 0.576-y) cos 0.987* + (0.030 - 1.152?) sin 0.987*] <= + .... 

Example I. 

If we put e = y — 0.1 and consider only the linear terms in «, these solu- 
tions become 

Mi = e-0.688< [0.1000 cos 0.987* - 0.0633 sin 0.987*], 
evi = - e-0.688< [0.0751 cos 0.987* + 0.0852 sin 0.987*], 

t«2 = e-0.688< [0.0100 cos 0.987* + 0.0701 sin 0.987*], 
&2 = e -0.688< [0.0830 cos 0.987* - 0.0085 sin 0.987*]. 



The values of the above quantities for various values of t are listed in Table II. 
The corresponding values for e u 3 and e v s have been computed from the center 
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of gravity equations. The coordinates (e «i, e v{), (e w 2 , « v 2 ) and (e m 3 , e %) de- 
note the x- and ^-displacements of mi, m% and m3, respectively, from the vertices 
of the equilateral triangle in configuration I. 



Table II 
= 0.1 7 = 0.1 



t 


e«i 


e»i 


ett 2 


e»2 


eM 3 


e#3 





+0.100 


-0.075 


+ 0.010 


+0.083 


-0.046 


-0.020 


0.1 


+0.087 


-0.078 


+0.016 


+0.076 


-0.044 


-0.015 


0.2 


+0.075 


-0.079 


+0.021 


+0.070 


-0.042 


-0.010 


0.3 


+0.063 


-0.079 


+0.024 


+0.063 


-0.040 


-0.006 


0.4 


+0.052 


-0.078 


+0.028 


+0 056 


-0 037 


-0.002 


0.5 


+0.041 


-0.076 


+0.030 


+0.049 


-0.034 


+0.001 


0.6 


+0.032 


-0.073 


+0.031 


+0.043 


-0 032 


+0.004 


0.7 


+0.023 


-0.069 


+0.032 


+0.036 


-0 029 


+0.006 


0.8 


+0.015 


-0.065 


+0.033 


+0.030 


-0.026 


+0.008 


0.9 


+0.008 


-0.061 


+0.033 


+0.025 


-0.023 


+0.010 


1 


+0.001 


-0.057 


+0.032 


+0.019 


-0.020 


+0.011 


1.2 


-0.009 


-0.047 


+0.030 


+0.010 


-0.014 


+0.013 


1.4 


-0.017 


-0.037 


+0.027 


+0.003 


-0.010 


+0.013 


1.6 


-0.021 


-0.028 


+0.023 


-0.003 


-0.005 


+0.013 


1.8 


-0.024 


-0.020 


+0.019 


-0.007 


-0.002 


+0.012 


2 


-0.025 


-0.012 


+0.015 


-0.010 


+0.001 


+0.011 


2.4 


-0.022 


-0.001 


+0.008 


-0.013 


+0.004 


+0.008 


2.8 


-0.017 


+0.006 


+0.002 


-0.012 


+0.005 


+0.005 


3.2 


-0.011 


+0.008 


-0.001 


-0.009 


+0.005 


+0.002 


3.6 


-0.006 


+0.009 


-0.003 


-0.006 


+0.004 


+0.0002 


4 


-0.0015 


+0.0072 


-0.0037 


-0.0033 


+0.0028 


-0.0009 


4.5 


+0.0015 


+0.0046 


-0.0032 


-0.0006 


+0.0011 


-0.0015 


5 


+0.0027 


+0.0021 


-0.0021 


+0.0008 


+0.0002 


-0.0014 


5.5 


+0.0026 


+0.0003 


-0.0011 


+0.0014 


-0.0004 


-0.00097 


6 


+0.0019 


+0.0006 


-0.0003 


+0.0013 


-0.0006 


-0.00052 


7 


+0.0004 


+0.0009 


+0.0004 


+0.0005 


-0.0004 


+0.00005 



The diagram of the orbit of Wi for Example I is shown in Figure 3. This 
orbit is, of course, with respect to the rotating axes. The arrow indicates the 
directions of motion. 
Example II. 

We have also computed the orbits when e = 0.1 and 7 = 1. The equations 
determining the displacements, in so far as the linear terms in e are concerned, are 

mi = e -0.688( (0.100 cos 0.987/ - 0.228 sin 0.9870, 

61 = - e-o.688( (0.231 cos 0.987/ + 0.060 sin 0.987/), 
e% = e-o.688< (0.100 cos 0.987/ + 0.111 sin 0.987/), 

62 = e -0.68St (0.135 cos 0.987/ - 0.112 sin 0.987/). 



On substituting in these equations the various values of t as in Table II and mak- 
ing use of the center of gravity equations to determine e u 3 and e v^, the following 
displacements are found : 
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Fig. 3. 
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Fig. 4. 
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0.1 



= 1 



t 


««i 


e»i 


««2 


ei>2 


eu 3 


ev 3 





+0.100 


-0.231 


+0.100 


+0.135 


-0.100 


+0.011 


0.1 


+0.072 


-0.221 


+0.103 


+0.115 


-0.091 


+0.019 


0.2 


+0.046 


-0.208 


+0.105 


+0.096 


-0.081 


+0.026 


0.3 


+0.024 


-0.195 


+0.104 


+0.079 


-0.072 


+0.031 


0.4 


+0.004 


-0.179 


+0.103 


+0.062 


-0.063 


+0.034 


0.5 


-0.014 


-0.164 


+0.100 


+0.047 


-0.054 


+0.038 


0.6 


-0.029 


-0.149 


+0.096 


+0.033 


-0.046 


+0.040 


0.7 


-0.042 


-0.133 


+0.092 


+0.020 


-0.038 


+0.041 


0.8 


-0.053 


-0.119 


+0.086 


+0.008 


-0.031 


+0.042 


0.9 


-0.061 


-0.104 


+0.080 


-0.001 


-0.024 


+0.042 


1.0 


-0.068 


-0.090 


+0.074 


-0.010 


-0.018 


+0.041 


1.2 


-0.076 


-0.063 


+0.062 


-0.023 


-0.007 


+0.039 


1.4 


-0.078 


-0.039 


+0.049 


-0.032 


+0.002 


+0.035 


1.6 


-0.076 


-0.019 


+0.037 


-0.038 


+0.009 


+0.030 


1.8 


-0.071 


-0.003 


+0.026 


-0.040 


+0.013 


+0.025 


2 


-0.063 


+0.009 


+0.016 


-0.040 


+0.016 


+0.020 


2.4 


-0.044 


+0.024 


+0.001 


-0.034 


+0.017 


+0.011 


2.8 


-0.026 


+0.028 


-0.008 


-0.024 


+0.015 


+0.003 


3.2 


-0.011 


+0.026 


-0.011 


-0.015 


+0.011 


-0.001 


3.6 


-0.0001 


+0.0198 


-0.0114 


-0.0066 


+0.0069 


-0.0039 


4 


+0.0061 


+0.0130 


-0.0095 


-0.0008 


+0.0033 


-0.0047 


4.5 


+0.0087 


+0.0054 


-0.0061 


+0.0032 


+0.0001 


-0.0041 


5 


+0.0078 


+0.0003 


-0.0028 


+0.0045 


-0.0015 


-0.0028 


5.5 


+0.0054 


-0.0024 


-0.0004 


+0.0039 


-0.0019 


-0.0014 


6 


+0.0028 


-0.0031 


+0.0009 


+0.0027 


-0.0017 


-0.0004 


7 


-0.0004 


-0.0018 


+0.0012 


+0.0004 


-0.0005 


-0.0005 



The diagram of the orbits of the three bodies in this example is found in Fig- 
ure 4. The length of the side of the equilateral triangle is taken as the unit 
length. 

In conclusion, the author wishes to express his thanks to his former colleague, 
F. M. Wood, M.A., B.Sc, for verifying the algebraic and numerical computa- 
tions in this paper. 

The University op British Columbia, 
Vancouver, Canada. 



ERRATA, VOLUME 23 
Page 51. J. F. RiTT, Prime and composite folynominals. 
Page 63. Line 7, for ^~ 1 (v) read *p~ l (w); 

line 25, for ¥> -1 ( u ) read ip^iw). 



